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Quiescent X-ray variability from the neutron star transient 
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ABSTRACT 

A number of studies have revealed variability from neutron star low-mass X-ray bina- 
ries during quiescence. Such variability is not well characterised, or understood, but 
may be a common property that has been missed due to lack of multiple observations. 
One such source where variability has been observed is Aql X-l. Here, we analyse 14 
Chandra and XMM-Newton observations of Aql X-l in quiescence, covering a period 
of approximately 2 years. There is clear variability between the epochs, with the most 
striking feature being a flare-like increase in the flux by a factor of 5. Spectral fitting 
is inconclusive as to whether the power-law and/or thermal component is variable. 
We suggest that the variability and flare-like behaviour during quiescence is due to 
accretion at low rates which might reach the neutron star surface. 
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1 INTRODUCTION 

Transient neutron star low-mass X-ray binaries spend the 
majority of their lifetime in a quiescent state and only a 
small amount of the time in outburst, accreting at a sig- 
nificantly higher rate (~ 0.1 — — 1.0L e dd). During quies- 
cence the X-ray spectrum is typically characterised by a 
thermal and/or power-law component. The thermal com- 
ponent is most frequently interpreted as emission from the 
neutron star surface. In the deep crustal heating scenario 
l|Brown. Bildsten fc Rutl edgc 199cj), pycnonuclear reactions 
occur in the inner crust during outburst due to compres- 
sion of the crust by accreted material. Energy deposited by 
these reactions heats the neutron star core on a timescale of 
10 4 — 10 5 yr, where it reaches a steady state luminosity set 
by the time-averaged mass accretion rate. The neutron star 
should therefore have a minimum thermal luminosity due to 
these processes. This thermal emission is fit well by neutron 
star atmosphere models, which could open the possibility 
of m easuring the neutron star radius (e.g., iRutledge et alj 
1 19991 ). 

The origin of the power-law component, however, re- 
mains poorly understood. A number of emission mechanisms 
have been suggested involving on-going low levels of accre- 
tion, but whether the gas impacts the neutron star surface 
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(e.g.. iMenou fc McC lintock 200 ll). is stopped at the magne- 
tospheric radius, or whether there is a shock between a ra- 
dio pulsar relativistic wind and ma tter transferred from the 
companion star (see the review bv lCampana et aI1ll998l ) is 
unclear. Moreover, on-going accretion onto the neutron star 
surface could potenti ally produce a thermal-like spectrum 
IjZampieri et al.ll 199a ). 

Variability has been observed during quiescence in 
a number of sources, with variability seen over a wide 



years ( 


Campana et al. 1997. 2004: Rutledge et al 


200C. 


2001b, 


2002: ICackett et al.l 


2005. 20101: iMuno et al.1 


2007; 


Fridriksson et al. 2010. 2011 


). In the majority of cases, the 



variability has either been attributed to the power-law com- 
ponent or its origin remains unclear. However, our recent ob- 
servations of Cen X-4 showed variability as large as a factor 
of 4.4 which could only be fit by v ariability in both the ther- 
mal and power-law components l|Cackett et al.l [2010). An- 
other source showing parti cularly interesting quiescent v ari- 
ability is XTE J1701-462 l|Fridriksson et alj|2010l , l201ll ). In 
this source, an overall cooling of the thermal component has 
been observed after a long outburst lasting ~2 years. How- 
ever, on two occasions a sudden short-term flare has been 
observed, with the brightest one rising to a flux about 20 
times the normal quiescent level. Such flares are presumably 
due to sporadic increases in accretion rate during quiescence. 

The focus of this current work is Aql X-l, one of 
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Figure 1. RXTE/ASM 1-day averaged lightcurve of Aql X-l 
around the times of the Chandra and XMM-Newton quiescent 
observations. Significant detections (3cr) are marked as red open 
circles, whereas non-detections are black filled circles. The blue 
stars mark the times of the quiescent observations with Chandra 
and blue diamonds mark the times of XMM-Newton observations. 

the sources where variability has been observed both on 
short (hundreds of se conds) and long (months) timescales 
ijRutledee et alJl2002h . How ever, the nature of thi s variabil- 
ity has been debated. While IRutledge et all (|2002l ) conclude 
that the thermal componen t is va riable, a subsequent anal- 
ysis by ICampana fc Stellal (2003) showed that correlated 
changes in the power-law index and column density can also 
describe the data. Having recently found large amplitude 
quiescent variability in Cen X-4, we also searched for further 
variability in archival observatio ns of Aql X-l. In add ition 
to the 4 observations ana lysed bv lRutledee et all (|2002h and 
ICampana fc Stellal l|2003l ). we analyse a further 10 quiescent 
observations performed about a year later. In section [2] we 
detail the data reduction, section [3] describes the spectral 
analysis, and in section |4] we discuss our findings. 



2 DATA REDUCTION 

Here we analyse a total of 14 observations of Aql X-l in a 
quiescent state (11 with Chandra, 3 with XMM-Newton), all 
the available archival Chandra and XMM-Newton observa- 
tions of Aql X-l in quiescence. In Tabled] we give details of 
these 14 observations. Figure [T] shows the one-day averaged 
lightcurve from the RXTE All-sky monitor (ASM), with the 
times of these observations marked. Two outbursts occur 
between the first four Chandra observations and the last 10 
Chandra/ XMM-Newton observations. All observations anal- 
ysed here are clearly at times when Aql X-l is undetected 
by the RXTE/ASM. The 1-day average detection thresh- 
old (2 - 10 keV) fo r the ASM is approximately 10 mCrab 
|Levine et al.lll996l ). 

2.1 Chandra data reduction 

We reduced the Chandra data using the CIAO software 
(v4.2) and the most recent calibration database (CALDB 
v4.3.1). All observations were performed using the ACIS- 
S detector, with a 1/8 sub-array to give a frame time of 



0.44s and the source placed off-axis where the PSF is sig- 
nificantly larger. Bot h these choices mitigate pile-up (see 
IRutledge etall 2002), even for the brightest observation 
here. In all cases, we followed the standard data reduction 
thread^], using the chandra_repro script to reprocess the 
data with the latest calibration. The source spectrum was 
extracted using a circular region of radius 10 pixels, and the 
background spectrum was extracted from a source-free an- 
nulus with inner radius 20 pixels and outer radius 60 pixels. 
The response matrix (rmf) and ancillary response file (arf) 
were generated with the mkacisrmf and mkarf scripts. 

2.2 XMM-Newton data reduction 

All three XMM-Newton observations analysed here were per- 
formed with the PN detector operated in timing mode. As 
the source was in quiescence during these observations, the 
timing mode data is not of a high quality and the source is 
only detectable when a restricted energy range is used. We 
therefore analyse only the data for the two MOS detectors, 
which were both operated in full frame mode with the thin 
filter for each observation. These MOS data were reduced 
using the XMMSAS software (vlO.0.0), with the latest 'cur- 
rent calibration files', producing calibrated event lists from 
the Observation Data Files using emproc. 

There was no significant background flaring in the first 
(0112440301) or second (0112440401) XMM-Newton obser- 
vations analysed here. However, during the third observa- 
tion (0112440101), there was a bright background flare at 
the end of the observation. We searched for background 
flares by filtering the lightcurve from the whole detector 
for single events (pattern zero) with energies above 10 keV. 
We excluded times when the count rate from this filtered 
lightcurve was greater than 2 counts s . This reduced the 
exposure by approximately 0.5 ks. 

Using the evselect tool, the source spectra were ex- 
tracted from a circular region of radius 32", where as the 
background spectra were extracted from an annulus with in- 
ner radius 60"and outer radius 250". We selected for events 
with patterns 0-12 only. The rmf and arf were generated 
with the rmf gen and arf gen tools. 



3 SPECTRAL ANALYSIS 

The spectra we re modelled us ing the XSPEC (vl2) spectral 
fitting package |Arnaud| [l996). All spectra were grouped to 
a minimum of 20 counts per bin in the 0.5 - 10 keV energy 
range, and the models were only fit over this energy range. 
All uncertainties are quoted at the lo confidence level. 

An initial look at the raw count rates from the eleven 
Chandra observations (see Table [TJ demonstrates there is 
variability between the observations. For instance, CX07 
has a count rate about 4 times higher than CXO10. Initially 
we fit all the spectra individually with an absorbed neu- 
tron star atmosphere model, using the phabs model for pho- 
toele ctric absorption and the nsatmos model l|Heinke et al.l 
2006) for the neutron star atmosphere. From these fits, it is 
clear that half of the spectra are fit well {xt < L2) by just 

1 see http://cxc.harvard.edu/ciao/threads/ 



© 0000 RAS, MNRAS 000, 000-000 



Quiescent variability from Aql X-l 3 

Table 1. Quiescent observations of Aql X-l 



Observation 


OLetrL QctTiC 


1V1JU 


Mission/ Instrument 


V ) IJbLLJ 


Exp. time 


Net count rate 




dd/mm/yy 


mid observation 






( kscc) 


(c s J 


CXOl 


28/11/00 


51876.5 


Chandra/ ACIS-S 


708 


6.6 


0.182 ± 0.005 


CX02 


19/02/01 


51959.5 


Chandra/ACIS-S 


709 


7.8 


0.097 ± 0.004 


CX03 


23/03/01 


51991.9 


Chandra/ACIS-S 


710 


7.4 


0.129 ± 0.004 


CX04 


20/04/01 


52019.6 


Chandra/ACIS-S 


711 


9.2 


0.124 ±0.004 


CX05 


04/05/02 


52399.0 


Chandra/ACIS-S 


3484 


6.5 


0.166 ± 0.005 


CX06 


20/05/02 


52414.4 


Chandra/ACIS-S 


3485 


7.0 


0.187 ±0.005 


CX07 


11/06/02 


52436.2 


Chandra/ACIS-S 


3486 


6.5 


0.349 ± 0.007 


CX08 


05/07/02 


52460.7 


Chandra/ACIS-S 


3487 


5.9 


0.100 ±0.004 


CX09 


22/07/02 


52477.9 


Chandra/ACIS-S 


3488 


6.5 


0.091 ± 0.004 


CXO10 


18/08/02 


52504.4 


Chandra/ACIS-S 


3489 


7.1 


0.083 ± 0.003 


CXOll 


03/09/02 


52520.7 


Chandra/ACIS-S 


3490 


6.9 


0.107 ±0.004 


XMM1 


15/10/02 


52562.1 


XMM-Newton/MOS 


0112440301 


7.1 


0.051 ±0.003 (MOS1) 
0.055 ±0.003 (MOS2) 


XMM2 


17/10/02 


52564.2 


XMM-Newton/MOS 


0112440401 


13.4 


0.047 ± 0.002 (MOS1) 
0.047 ± 0.002 (MOS2) 


XMM3 


27/10/02 


52574.1 


XMM-Newton/MOS 


0112440101 


2.7 


0.057 ±0.005 (MOS1) 
0.057 ±0.005 (MOS2) 



an absorbed neutron star atmosphere model, but that the 
other half require an additional power-law component. The 
observations that require a power-law component (Xv > 1-5 
without it) are CX03, CX04, CX05, CX06, CX07, CXOll 
and XMM2. Note that in most of those fits the power-law 
index is poorly constrained, and all spectra are consistent 
with having the same index at the la level. In these fits 
where the power-law is required, the index ranges from ap- 
proximately -1 to 2, but the uncertainties are such that they 
are all consistent with T = 1. 

In order to determine which spectral parameters are 
variable between the observations, we fit all the spectra 
jointly. Here, we fit an absorbed neutron star atmosphere 
plus power-law model. We include the power-law component 
for all observations, allowing the normalisation to be a free 
parameter. Given that in the individual fits the power-law 
index was consistent with being constant, we tied the power- 
law index between all observations. From the individual fits, 
we also find that the column density is consistent between 
all the observations, and therefore we tie this parameter. 
We also assume a canonical neutron star mass and radius 
(1.4 Mq, 10 km), and that the entire neutron star surface 
is emitting (these parameters were all fixed). We assume a 
distance of 5 kpc throughout, as also adopted in previous 
work (e.g. iRutledee et"al]|2002h. For a discussion of the dis- 
tance to Aql X-l see lRutledge et al.l ^OOlal ) who find the 
distance is between 4 and 6.5 kpc. 

In the first instance, we tie the neutron star atmosphere 
temperature between all observations (the spectral parame- 
ters are given in Table [2] and all uncertainties quoted there, 
and throughout the paper, are at the la confidence level). 
This allowed us to investigate whether variability in the 
power-law normalisation alone can fit the data. This pro- 
vides an adequate fit to the data (x„ = 1.10, v — 533, null 
hypothesis probability = 0.057), with an effective tempera- 
ture (for an observer at infinity) of fcT c g = 108.9 ± 0.8 eV, 
power-law index, F — 2.73 ± 0.06 and power-law normali- 
sation ranging from a non-detection in XMM2 to approxi- 
mately 10" 3 photons keV -1 cm~ 2 s _1 at 1 keV in CX07. 
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Figure 2. Variability of the unabsorbed 0.5 - 10 keV power-law 
flux when this is the only parameter allowed to vary between 
observations (the effective temperature and power-law index are 
tied between the observations). 



The variability of the unabsorbed 0.5 - 10 keV power-law 
flux is shown in Figure [2] for this model. Note that the col- 
umn density we determine is consistent with values deter- 
mined from previous work, as well as optica l photometry and 
21 cm emission (see iRutledge et al. 2002, for a discussion 
of the reddening towards Aql X-l). 

While the model with only the power-law normalisation 
variable between observations fits the data well, examining 
the residuals of the brightest observation (CX07) it seems 
that the steep power-law is attempting to mimic a soft ther- 
mal component, and the residuals show an upturn above 
4 keV (see Figure [3} • Note that while this is not statisti- 
cally significant, it would appear to suggest the spectrum 
should be dominated by a thermal component at soft ener- 
gies rather than the power-law. We therefore also fit the data 
allowing both the power-law normalisation and the temper- 
ature of the neutron star atmosphere to vary. The results of 
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Table 2. Spectral fitting parameters. The spectra were fit jointly, with the column density and power-law index tied between all 
observations in all cases. The neutron star mass and radius were fixed at canonical values (1.4 Mq, 10 km). The power-law normalisation 
is in units of photons keV -1 cm~ 2 s _1 and is defined at 1 keV. The luminosity, L, is evaluated over the 0.5 — 10 keV range and is for 
D = 5 kpc. The uncertainty in the luminosity is generally dominated by the uncertainty in D. 



Obs. 


N H 

(10 21 cm" 2 ) 


(eV) 


r Power law norm. 0.5 — 10 keV flux 

(10~ 5 ) (10~ 13 erg cm -2 s -1 ) 


L 

(10 33 erg s _1 ) 


Thermal 
fraction 


xl (dof) 


Power-law normalisation variable, temperature tied 



CXOl 


4.3 ± 0.1 


108.9 ± 0.8 


2.73 ± 0.06 


34.2 ± 3.2 


8.0 ± 0.6 


5.6 ± 2.7 


0.41 ± 0.03 


CX02 








7.0 ± 2.0 


3.8 ± 0.4 


2.9 ± 1.4 


0.77 ± 0.09 


CX03 








18.3 ± 2.3 


5.5 ± 0.5 


4.0 ± 2.0 


0.56 ± 0.05 


CX04 








18.6 ± 2.1 


5.6 ± 0.5 


4.1 ± 2.0 


0.56 ± 0.05 


CX05 








34.5 ± 3.1 


8.0 ± 0.7 


5.6 ± 2.7 


0.40 ± 0.04 


CX06 








45.3 ± 3.3 


9.7 ± 0.6 


6.7 ± 3.2 


0.34 ± 0.02 


CX07 








107.8 ± 5.6 


19.4 ± 0.9 


12.7 ± 6.1 


0.18 ± 0.01 


CX08 








10.7 ± 2.3 


4.4 ± 0.6 


3.3 ± 1.6 


0.68 ± 0.10 


CX09 








7.8 ± 1.9 


3.9 ± 0.5 


3.0 ± 1.5 


0.75 ± 0.10 


CXO10 








3.5 ± 1.9 


3.3 ± 0.5 


2.6 ± 1.3 


0.87 ± 0.11 


CXOll 








13.9 ± 2.1 


4.9 ± 0.6 


3.6 ± 1.8 


0.63 ± 0.08 


XMM1 








3.3 ± 1.8 


3.2 ± 0.4 


2.6 ± 1.3 


0.87 ± 0.11 


XMM2 








5+ 14 
U -0.5 


2.8 ± 0.4 


2.3 ± 1.2 


97+ 03 
u ' a ' -0.14 


XMM3 








3.6 ± 2.5 


3.3 ± 0.5 


2.6 ± 1.3 


0.86 ± 0.13 



1.10 (533) 



Temperature and power-law normalisation variable 



CXOl 


3.8 ± 0.1 


126 


4 ± 


.9 


0.80 ± 0.25 


0.9 ± 0. 


3 


8.0 ± 0.8 


4.9 ± 2.4 


0.89 ± 0.09 


CX02 




109. 


3 ± 1 


.3 




1.2 ± 0. 


6 


5.2 ±2.1 


3.1 ± 1.9 


76+ ' 24 
° ' -0.31 


CX03 




116 


1 ± 1 


.0 




1.4 ± 0. 


4 


6.8 ± 1.2 


3.9 ± 2.0 


0.78 ± 0.14 


CX04 




114 


.6 ± 1 


.0 




1.9 ± 


7 


7.5 ± 1.2 


4.0 ± 2.0 


0.71 ± 0.12 


CX05 




124 


2 ± 1 


1 




1.7 ± 0. 


6 


9.0 ± 1.3 


5.1 ± 2.6 


0.80 ± 0.12 


CX06 




125. 


2 ± 1 


1 




3.5 ± 1. 


3 


12.6 ± 2.0 


6.4 ± 3.2 


0.67 ± 0.11 


CX07 




142 


4 ± 1 


.3 




7.2 ± 2. 


6 


25.1 ± 4.1 


11.8 ± 6.0 


0.63 ± 0.11 


CX08 




112. 


2 ± 1 


.4 




1.0 ± 0. 


6 


5.3 ± 2.1 


3.2 ± 2.0 


n o,+0.19 


CX09 




110. 


1 ± 1 


.2 




0.9 ± 


3 


4.9 ± 1.0 


3.0 ± 1.6 


0.80 ± 0.17 


CXO10 




107 


1 ± 1 


.4 




1.2 ± 0. 


5 


5.0 ± 1.6 


2.8 ± 1.6 


0.74 ± 0.24 


CXOll 




113 


4 ± 1 


1 




1.2 ± 0. 


4 


6.0 ± 0.9 


3.5 ± 1.8 


0.79 ± 0.12 


XMM1 




107 


9 ± 1 


.0 




0.6 ± 0. 


2 


3.9 ± 0.8 


2.5 ± 1.3 


°-MtS:ii 


XMM2 




106 


2 ± 


.9 




0.5 ± 0. 


2 


3.5 ± 0.7 


2.3 ± 1.2 


0.88+-- 


XMM3 




109 


9 ± 1 


.1 




o.o+°- ! 


i 


3.1 ± 1.7 


2.4 ± 1.7 


l-00_o.55 



those fits are also given in Table [2] A good fit, with a lower 
X 2 , is achieved (xl = 0.88, v — 520). In this case a much flat- 
ter power-law index is found, Y — 0.80 ± 0.25. A power-law 
consistent with a slope of 1 - 2 is quite typical for quiescent 
neutron stars. With this model we find significant variability 
in the effective temperature of the neutron star atmosphere 
component, ranging from approximately kT^§ = 106 — 142 
eV (with typical la uncertainties of 1 eV), in addition to 
variability in the power-law normalisation. In Figure [4] we 
show the 0.5 - 10 keV flux, effective temperature and un- 
absorbed 0.5 - 10 keV power-law flux from this model. For 
comparison, we show the spectral fit to CX07 for both mod- 
els in Figure [3] note that the residuals are much flatter at 
both high and low energy in the fit where both the temper- 
ature and the power-law normalisation are allowed to vary. 

We also investigated a fit with both the power-law index 
and normalisation tied between all observations. It provides 
a significantly worse fit (xl = 1.16, v — 533, null hypothesis 
probability = 6.6 x 10~ 3 ), and cannot match the spectrum of 
the brightest observations above 3 keV, underpredicting the 
flux at 5 keV by a factor of ~ 4. A s discussed in the Intro- 
duction, ICampana fc Stellal l|2005tl described the first four 
Chandra observations by correlated changes in the power- 
law index and column density. Note that when fitting the 
spectra individually, the Nu values are all consistent with 
remaining unchanged, and as mentioned above, so are the 
power-law indices. Moreover, given that a good fit can be 



achieved without the power-law index or JVh being variable, 
we do not test this model further. 

The three XMM-Newton observations have some of the 
lowest 0.5 - 10 keV fluxes of all the observations. While this 
change may be real, we caution that part of the difference 
in flux with the closest Chandra observations could poten- 
tially be due to cross-calibration differences between the two 
missions. Differences of 10% between Chandra and XMM- 
Newton soft X-ray fluxes have been seen in cross-calibration 
studies (e.g. iTsuiimoto et al.ll2oilT ). 

We also investigated the effect of having fit all the data 
simultaneously with multiple parameters tied between the 
spectra. With such a method, the spectrum with the high- 
est signal-to-noise ratio can potentially dominate the fit, and 
skew the parameters. However, we find that this is not the 
case here. Removing the best spectrum, CX07, from the fit 
with both temperature and power-law normalisation vari- 
able, we find that both tied parameters (the power-law in- 
dex and Ah) remain consistent with their former values, 
with T = 0.79 ± 0.35 and A H = (3.8 ± 0.1) x 10 21 cm -2 . 

As an aside, we note that we searched for any obvi- 
ous variability or flaring within the individual lightcurves of 
each observ ation. There is no dr amatic variability present 
(though see iRutledge et al.l |2002| . for power density spec- 
tra of the first four Chandra observations), and we do not 
investigate this further here. 
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Figure 3. The brightest Chandra spectrum of Aql X-l, CX07. 
Panel (a) displays the model from jointly fitting all spectra with 
only the power-law normalisation variable. Panel (b) shows the 
model from jointly fitting all spectra with both the effective tem- 
perature of the neutron star atmosphere and the power-law nor- 
malisation variable. In both cases the dashed red line marks the 
neutron star atmosphere component, and the dotted blue line 
marks the power-law component. 



4 DISCUSSION 

We have analysed 14 observations of Aql X-l in quiescence, 
with the observations spanning a period of approximately 
2 years. Variability in the first four observations has been 
studied previously (|Rutledge et al.ll2002l ; ICampana fc Stellal 
2003). Here, we find variability between the different obser- 
vations, with a particularly striking flare showing a change 
in 0.5 - 10 keV flux by a factor of ~5 between the peak of 
the flare (0XO7) and the next observation (0XO8). If the 
three observations around MJD 52400 are from the same 
event, then the flare lasted approximately 60 days, with ob- 
servations CX05 and CX06 seemingly catching the rise of 
the flare. However, with the sparse sampling we cannot be 
sure that this was a single flare event or not. The ratio of 
the maximum to minimum observed 0.5 - 10 keV flux is 7 - 
8 (depending on the spectral model assumed). 

We fitted the spectra with an absorbed neutron star 
atmosphere plus power-law model. In the first instance, 
the only parameter allowed to vary between the observa- 
tions was the power-law normalisation, i.e. we assume that 
the temperature of the neutron star surface remains the 
same at all epochs. Such a model fits the data adequately 



— 30 




52600 



Figure 4. Quiescent lightcurve of Aql X-l for the model where 
temperature and power-law normalisation are allowed to vary. 
The top panel shows the 0.5 - 10 keV flux, the middle panel the 
effective temperature (for an observer at infinity), and the bottom 
panel the unabsorbed 0.5 - 10 keV power-law flux. 



(xl = 1-10), and requires a power-law with a steep slope 
(r = 2.7). Alternatively, we also fit the data allowing both 
the effective temperature and power-law normalisation vary, 
which, of course, also fits well (\l = 0.88). In this case, the 
power-law index is much flatter (F = 0.8), and the thermal 
component dominates the spectrum at all epochs (with ther- 
mal fractions typically 0.7 - 0.9). Given that both models 
provide a good fit to the data, we are unable to conclusively 
determine which spectral component is driving the variabil- 
ity. 

Comparing the times of the quiescent observations with 
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Figure 5. The RXTE/ASM 1-day averaged lightcurve of Aql X- 
1 around the times of the three brightest quiescent observations. 
There is no apparent flaring detected in the ASM lightcurve. The 
markers are as in Figure [T] 

the RXTE/ASM lightcurve (Fig [[J, the last significant (3<r) 
detection at the end of the outburst (before the flare) with 
the RXTE/ASM was on MJD 52355. A pointed RXTE/PCA 
observation detects Aql X-l on MJD 52366 with a 2 - 10 
keV flux of 8.5I3 3 x 10 -12 erg s _1 cm -2 , however, a sub- 
sequent (and final) pointed observation two days later does 
not provide a clear detection, and gives an upper limit to 
the 2-10 keV flux of approximately 5.5 x 10~ 12 erg s" 1 
cm" 2 . Therefore, the first of the three bright observations 
(CX05, 6 and 7) occurs about 33 days after the end of a 
full outburst (taking the last detection with the PCA as the 
approximate end of the outburst). 

Figure [S] shows the RXTE/ASM 1-day averaged 
lightcurve around the times of the three brightest quiescent 
observations. Clearly, there is no apparent flaring detected 
by the ASM. There is only one data point at a significance 
greater than 3<r, and that occurs approximately 5 days be- 
fore CX05 (the observation on MJD 52399). Given the sur- 
rounding non-detections, this point is likely a statistical fluc- 
tuation and not a real detection. Even so, Aql X-l was not 
accreting at close to outburst levels near the quiescent flare 
we observed. 

One source of variability in quiescence is thermal re- 
laxati on of the neutron sta r crust after the end of an out- 
burst ( Rutledge et al. 20 021), which has now been observed 
in four sources (e.g. ICackett et alj 120061 ; iFridriksson et ail 
l20ld ; iDegenaar et alJl201Ch . As CXQ5. 6 and 7 occur close 
to the end of an outburst, it may be possible that thermal re- 
laxation of the crust could lead to variability in the thermal 
component. Thus far, only an overall luminosity decrease 
over time has been observed. However, depending on the mi- 
crophysics of the crust and core, an increase can be achieved 
as hea t from pycnonuclear reactions diffuses to the top of th e 
crust (|Ushomirskv fc Rutledgell200ll ; iRutledge et~aii r2002). 
The largest fractional variability would come from sources 
with long recurrence times between outbursts and from neu- 
tron stars with enhanced levels of core neutrin o emission. 
But, as noted bv lUshomirskv fc Rutledge! (|200^ . the quies- 
cent luminosity of Aql X-l is inconsistent with rapid cool- 
ing I Brown et al.|[l99^ 1 and it has a short recurrence time 



(~ 1 yr), thus the observed quiescent variability in this 
source cannot be due to thermal relaxation of the crust. 
Furthermore, the observed flare is not just due to the ther- 
mal component - the thermal fraction is lowest during the 
flare and there is a significant and clearly variable power- 
law component. Note also that CXOl also occurs directly 
after the end of an outburst, and is hotter than the 3 sub- 
sequent observations. While it is tempting to speculate that 
the neutron star was hotter in CXOl than CXO 2-4 
b ecause of thermal relaxat i on of the crust, the arguments 
of lUshomirskv fc Rutledge] |200ll ) (discussed above), would 
seem to rule that out. 

While thermal relaxation of the crust does not appear to 
explain the flare, it may still be possible that the flare is as- 
sociated with the previous outburst. A double outburst was 
recently seen from the accretion-powered millisecond pul- 
sar ICR J00291+59 34 (|Patrundl2o"Tol ; lHartman et al.ll201ll ; 
iPapitto et ai]|201lt ). where two approximately 10 day long 
outbursts were separated by o nly 30 days in quiesc ence. 
Discussing these two outbursts, lHartman et al.l l |201ll ) sug- 
gest that they are connected - the first outburst is stopped 
by the propeller effect before the disc is completely de- 
pleted, leaving material in the disc for the second outburst. 
The first brighter observation in Aql X-l is also approxi- 
mately 30 days after the end of an outburst, however, un- 
like ICR J00291+5934 the increase in Aql X-l does not 
get close to outburst levels (Aql X-l is not detected by the 
RXTE/ASM near the brightest observations), thus it is not 
clear whether a similar mechanism applies here. 

In addition to Aql X-l, flaring (by which we mean short 
lasting, sharp rises in flux) during quiescence has recently 
been observed in the neutron star low-mas s X-ray bina- 
ries XTE J1701-46 2 (|Fridriksson et alJl201Cft and Cen X-4 
|Cackett et alj|2010h . In XTE J1701-462, an increase in X- 
ray flux by a factor of almost 3 (to a peak observed 0.5 - 
10 keV luminosity of 2.65 x 10 34 erg s _1 ) was observed ap- 
proximately 230 days after the end of an outburst on top of 
an overall decreasing X-ray flu x presumably associated with 
thermal relaxation of the crust (|Fridriksson et al . 2010). The 
increased flux lasted a maximum of approximately 120 days 
(the time between the quiescent observations on either side 
of the increased flux). Since then, further X-ray monitoring 
of this source during quiescence showed another X-ray flare 
about 1075 days into quiescence. This time, the flare was sig- 
nificantly brighter, with a peak observed increase of a factor 
of ~20 in flux (reaching a maximum lumin osity of ~ 1 x 10 35 
ergs -1 ) and lasting 10 - 20 days in total (|Fridriksso n et all 
120111) . Interestingly, like Aql X-l, in XTE J1701-462 the 
thermal component does increase during the flares, but the 
non-thermal component increases faster leading to a ther- 
mal fraction that is lower during the flare s than at other 
times in quiescence l|Fridriksson et al]|201ll ). 

In Cen X-4, quiescent observations spanning 15 
years showed subs tantial variability between epochs 
|Cackett et al.l l2010h . Additionally, examination of the 
lightcurves from individual observations revealed a short 
flare in one of the XMM-Newton observations. This flare 
lasted ~2 ks and showed an increase of approximately a 
factor of 4 in count rate at the peak. The natural explana- 
tion for this flaring activity in both these sources appears 
to be increased levels of ongoing accretion during quies- 
cence, which presumably is occurring in Aql X-l too. The 
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exact mechanism for this accretion is unclear (as discussed 
in the Introduction, the origin of the power-law component 
is uncertain), but the clear variability in the thermal com- 
ponent in Cen X-4, and the apparent connection between 
the power-law and thermal components (the thermal frac- 
tion remained roughly constant) may suggest that accreting 
material reaches the neutron star surface. 

Another class of objects that are particularly interest- 
ing are the very faint X-ray transients (VFXTs) whose peak 
2 

erg s 

missed by all-sky monitors. However, frequent monitoring of 
the Galactic Centre region by multiple X-ray missions (e.g. 
Chandra, XMM-Newton and Swift ), has reve a led a number 
of th e se sources (see, e.g. Muno et alj 120051: ISakano et ail 
2005 1; IWiinands et~ai1 l200rj ; iDegenaar fc Wiinandsl 120091. 
2010J). The frequent monitoring has revealed flare-like be- 



10 keV luminosity during outburst only reaches 10 34 36 
1 . With such low peak luminosities, these sources are 



haviour in-b etween normal quiescent and outbu rst levels in 
four sources i|Degenaar fe Wiinands|[2009l . |2010T | . The flares 
last typically 1-2 weeks, and the peak brightness can vary 
substantially, though one source (XMM J174457-2850.3) 
has shown a flare with a peak luminosity as low as seen 
here in Aql X-l. Several (though not all) of these VFXTs 
are confirmed as neutron star sources as they have shown 
type-I X-ray bursts. 

Variability in quiescence is not limited to just neu- 
tron stars. Several of the brighter quiescent stellar-mass 
black holes have also been observed to vary in several 
sources (V404 Cy g, 4U 1630-47, GX 3 3 9-4, A0620-00, and 



GRO J1655-40: iWagner et al.l Il994j; iParmar et~ai1 Il997l; 
Kong et alj |2000|. 120021; iHvnes et"ail |2004| ; iBradlev et alj 



200 



Miller- Jones et al.1 120081 ) . The best studied of these 
is V404 Cyg. This object has been observed to vary by 
as much as a factor of 20 in X-ray count rate during one 
short flare l|Hvnes et al. 2004|), and has been variable dur- 
ing the observation on eve ry occasion when observed by 
Chandra or XMM-Newto n l|Kong et all |2002| ; IHvnes et all 
120041 ; IBradlev et"aT1l2007l ). Also interesting is that the X- 
ray variability is correlated with optical variability, with 
the Ha line variability se emingly powered by X-ray irra- 
diation (|Hvnes et al.ll200"4h . Moreover, sensitive radio obser- 
vations of V404 Cyg in quiescence detect the source and 
even find a short flare, with a rise time of ~30 minutes 
(|Miller- Jones et al.ll200Sr ). The radio emission appears to be 
non-thermal and therefore may originate in a compact jet. 
Of course, in these objects there is no stellar surface, thus 
variability here cannot involve any stellar surface or stel- 
lar magnetic field. While it has been suggested that coro- 
nal emission from the rapidly rotating secondary star could 
be the s ource of the quiescent X-ray emission in black hole 
sources (|Bildsten fc Rutleded I2OO0I ). the X-ray luminosity 
of several sources exceeds the maxim um predicted by the 
coronal model (e.g. iKong et al.ll2002l ). 

Quiescent variability and flaring has also been seen in 
nearby quiescent supermassive black holes. For instance, the 
supermassive black hole at the centre of the Galaxy, Sgr A*, 
has shown significant fl aring behaviour at X-ray, IR and 



radio wavelengths (e.g. Baganoff et al 



l2003l : IZhao et al.1 120031 ; iPorquet et al 



2001; Gc nzel et al.1 



20031 . |200S| ). The 



brightest X-ray flare from Sgr A* showed an increase by 
a factor of 160 in 2- 10 keV luminosity w ith a duration 
shorter than one hour l|Porauet et al.ll2003l ). More recently, 



the supermassive black hole at the centre of Andromeda, 
M31*, has also shown X-ray varia bility and flaring activity 
l|Garcia et al.l l2010l : iLi et alll201Ch . Thus , while any mech- 
anism for such fla res is not clear (e.g. lYuan et al.l |2004| ; 
iMaitra et al.|[2009l . and references therein), it may suggest 
that flaring is a standard property of accretion at low rates, 
though the stellar surface and magnetic field of neutron stars 
may lead to a different inner accretion geometry. 

The standard model used to explain the transient be- 
haviour and outburst cycles in X-ray bin a ries is the disc 
instability model (e.g. King & Rittcr 1998; Hamcurv et al.l 
Il998l ; iMenou et al.l 12000 ; Du bus et al1l200ll ; lLasotal 120011 ). 
During quiescence, as long as the mass accretion rate is 
smaller than some critical value everywhere in the disc, 
accretion can continue (see for ex ample the discussion in 
iKuuIkers. in 't Zand fc Lasotall2009h . However, the disc in- 
stability model does not make a ny strong prediction s for 
variability during quiescence, and Dubus et al.1 l|200lf ) dis- 
cuss how flaring after the en d of an outburst such as that 
observed in GRO J0422+32 (|Callanan et al.lll995l ') are not 
understood. Similarly, it appears that the variability seen in 
Aql X-l and XTE J1701-462 can not readily be explained 
in the disc instability framework. 

Finally, the short 2 ks flare seen in the X-ray lightcurve 
of Cen X-4 may have a different origin (not related to the 
disc instability model). Such short 30 - 60 mi n flares are also 
seen i n the optical li g htcur ve of Cen X-4 (jShahbaz et al.l 
|2010| ). IShahbaz et~aH (|2010l ') suggest that the most likely 
model for those flares is blackbody radiation from an op- 
tically thin layer of recombining hydrogen (similar to the 
model for solar flares), and must occupy a small (0.3 per 
cent) area of the disc. Further monitoring of quiescent X- 
ray binaries is required to understand the prevalence of such 
variability and flares. 
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